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> ABSTRACT 

^ , Titanium carbide (TiC) nanocrystals were recently proposed as the carrier of the 

\ mysterious 21 |Um emission feature observed in post-asymptotic giant branch stars, based 

\ on their close spectral match and the presolar nature of meteoritic TiC nanograins 

QQ ' (which reveals their stellar ejecta origin). But we show in this Letter that the Kramers- 

■ Kronig dispersion relations, which relate the wavelength-integrated extinction cross 

y—^ • section to the total dust mass, would impose a lower bound on the TiC mass. This 

. Kramers-Kronig lower limit exceeds the maximum available TiC mass by a factor of at 

, least ~50, independent of the absolute value of the ultraviolet /visible absorptivity of 

\ nano TiC. The TiC model is therefore readily ruled out by the Kramers-Kronig physical 

Qh| principle. 

o : 

^ , Subject headings: circumstellar matter — dust, extinction — infrared: stars — stars: 

^ ; AGB and Post-AGB — stars: individual (HD 56126) 

> : 

' !• Introduction 

A prominent broad (with a FWHM of ~ 2 /im) emission feature at about 21 ^m^ was discovered 
by Kwok, Volk, & Hrivnak (1989) in an analysis of the Infrared Astronomical Satellite (IRAS) Low 
Resolution Spectrometer (LRS) spectra of four carbon-rich post-asymptotic giant branch (AGB) 
stars, and subsequently confirmed by both ground-based and airborne Kuiper Airborne Observatory 
and ISO observations (see Kwok, Volk, & Hrivnak 1999 and references therein). So far, this feature 
has been detected in twelve post- AGB stars (commonly termed as the "21 ^m sources"; Kwok et al. 
1999)^, with little shape variation found between different sources. These stars have quite uniform 



^ High resolution spectra obtained with the Short- Wavelength Spectrometer (SWS) instrument (with a resolution 
of A/ A A 2000) on board the Infrared Space Observatory show that this feature actually peaks at ~ 20.1 ^tm (Volk, 
Kwok, & Hrivnak 1999). But for historical reasons, this feature continues to be referred to as the "21 feature". 

weak 21 /im feature was recently reported for three planetary nebulae (PNs; Hony, Waters, & Tielens 2001; K. 
Volk 2003, in preparation). 
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properties: they are mostly metal-poor carbon-rich F and G supergiants with strong infrared (IR) 
excesses and overabundant s-process elements (see Kwok et al. 1999). 

The carrier of this mysterious 21 fim feature remains unidentified, although many candidate 
materials have been proposed including iron oxides Fe203 or FC3O4, hydrogenated amorphous 
carbon, hydrogenated fullerenes, hydrogenated nanodiamond, amides (thiourea or urea OC [NH2]2), 
SiS2, oxygen-bearing side groups in coal (see Kwok et al. 1999 and references therein), and more 
recently titanium carbide (TiC) nanoclusters (von Helden et al. 2000), SiC (Speck &; Hofmeister 
2003), and stochastically-heated silicon core-Si02 mantle nanograins (Smith & Witt 2002; Li & 
Draine 2002). 

The nano-TiC model seems attractive. While bulk TiC does not show any noticeable feature 
near 20.1 /xm (Henning & Mutschke 2000), laboratory absorption spectra of TiC nanocrystals con- 
taining 27—125 atoms exhibit a distinct feature at ~20.1/^m, closely resembling the astronomical 
21 iim emission feature both in peak position, width, and in spectral details (von Helden et al. 
2000). This model further gains its strength from the identification of presolar TiC grains (with 
radii ~ 100 A) in primitive meteorites as inclusions embedded in micrometer-sized presolar graphite 
grains (Bernatowicz et al. 1996). 

Since Ti is a rare element, it is important to know whether the amount of TiC dust required 
to account for the observed 21 /xm feature is a reasonable quantity. Since the optical properties 
of nano TiC in the ultraviolet (UV) /visible wavelength range are unknown, one has to rely on 
those of bulk TiC as a starting point. Such an attempt has recently been made for the post-AGB 
star HD 56126 by Hony et al. (2003). They found that to explain the observed 21 fim feature, the 
TiC model requires a much higher UV/ visible absorptivity than that of bulk TiC (by a factor of 
>20) . However, we will demonstrate in this Letter how the TiC absorption integrated over a finite 
wavelength range places a lower limit on the TiC dust mass through the Kramers-Kronig dispersion 
relations. For HD 56126, this Kramers-Kronig lower limit exceeds the maximum available TiC mass 
by a factor of >50, independent of the absolute value of the UV/visible absorptivity. Therefore, 
increasing the UV/visible absorptivity is unable to avoid the TiC abundance problem. 

2. HD 56126: A Test Case 

HD 56126, a bright (visual magnitude ~ 8.3) post-AGB star with a spectral type of F0-5I, is 
one of the four 21 yum sources originally discovered by Kwok et al. (1989). Mid-IR imaging of this 
object at 11.9 /xm shows that its circumstellar dust is confined to an area of 1.2"-2.6" from the 
star (Hony et al. 2003). Detailed modeling of its dust IR spectral energy distribution suggested a 
dn{r)/dr ~ 1/r dust spatial distribution at rmin ^f< Vmax, where r is the distance from the star, 
rmin and rmax are respectively the inner and outer edge of the HD 56126 dusty envelope (Hony et 
al. 2003). If TiC nanograins are indeed present in HD 56126 and follow the distribution of the bulk 



3 



(hydrogenated) amorphous carbon dust, the total power absorbed by the TiC dust would be^ 

Elt (TiC) = m^j^^l ^"^^— Z;-^ r " ^) ^A^^ , (1) 

^ ^max ~ Vin ^912 A 

where mJ^*(j is the total mass of the TiC nanoparticle component, Kabs{a, A) is the mass absorption 
coefficient (cm^ g~^) for nano TiC grains of size a at wavelength A, ry^ is the stellar radius, is the 
flux per unit wavelength (ergs~^ cm~^ /um"-'^) at the top of the illuminating star's atmosphere. In 
the nanometer size domain, TiC grains are in the Rayleigh limit and hence /«abs(^) is independent 
of size a in the wavelength range where the stellar radiation peaks. Following Hony et al. (2003), 
we will adopt a distance of d ~ 2.4 kpc to the star (and therefore rmin ~ 4.3 x lO"*^^ cm, rmax ~ 
9.3 X 10^^ cm), a stellar radius of ~ 49.2 (r© is the solar radius), a stellar luminosity of 
6054^0 {Lq is the solar luminosity), and approximate the HD 56126 stellar radiation by the 
Kurucz (1979) model atmospheric spectrum with T^s = 7250 K and log^ = 1.0. 

The total power emitted in the 21 fim feature of HD 56126 was estimated to be El'^ (21 /xm) « 
1.0 X lO^^ergs"^ (Hony et al. 2003). A lower limit on the total mass of the TiC dust required to 
account for the observed 21 /xm feature can be obtained by assuming that all the energy absorbed 
by the TiC dust would be emitted solely in this feature: 

2 ^2 _ „2 foo 

m^l^ = El^ (21/.m) / / „ K,bs(A, a) F^dA . (2) 

^i, in [ruiax/ruiin) 7912 A 

Apparently, the required TiC dust mass m™g is sensitive to the material's UV/visible absorptivity. 

On the other hand, an upper limit on the total mass of the TiC dust which could be present in 
the circumstellar envelope around HD 56126 can be estimated from its photospheric Ti abundance 
Ti/H 1.3 X 10~^ (van Winckel &; Reyniers 2000) and the envelope mass (including both gas and 
dust) 0.16 < menv < 0.44 m© (Hony et al. 2003): m^^ pa 2.5 x 10"'^ m©. 



3. Ultraviolet/ Visible Absorption Properties of TiC Nanograins 

Since the UV/visiblc absorption properties of nano TiC are unknown, we will first represent 
them by those of bulk TiC. Using Mic theory and the optical constants measured for bulk TiC 
samples (Koide et al. 1993), wc calculate the UV/visible absorption efficiency Qabs of spherical 
TiC nanograins. For illustration, we plot in Figure 1 the mass absorption coefficients (/«abs = 
SQahs/ [4apTic] where pxiC ~ 4.92 gcm~^ is the TiC mass density) of TiC spheres of radii a=10A 
calculated from the optical constants of bulk TiC. Using these Kabs values, we derive a lower limit 



^The circumstellar envelope around HD 56126 is assumed to be optically thin. If there actually exists appreciable 
extinction (e.g., Hony et al. [2003] found an average visual extinction of ^1.1 mag), the conclusion of this Letter 
would be strengthened since when exposed to an attenuated stellar radiation field, the TiC model would require more 
TiC dust to account for the same 21 fim feature strength. 
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Fig. 1. — Mass absorption coefficients Kabs(-^) of TiC nanograins of radii a = 10 A calculated (1) from tlie optical 
constants of bulk TiC (solid line), (2) from the Hony et al. (2003) hypothetical formula with K^bs = 5.6 x 10® cm^ 
(see Eq.[3]; dashed line), or (3) from a hypothetical Drude profile with k^j^ = 9.6 x lO^cm^g"^ (see §4; dot- 
dashed line). The IR part is the laboratory measured 21 fim spectrum of TiC nanocrystals (von Helden et al. 2000) 
approximated by a Lorentz profile (Chigai et al. 2003). It is apparent that the Hony et al. (2003) and Drude 
hypothetical absorption spectra greatly enhance the "ability" of TiC nanograins in absorbing the HD 56126 starlight. 



of tti^Iq ~ 9.3 X 10~^ m0 on the TiC dust mass from Eq.(2). Apparently, if nano TiC grains do 
not have a much higher UV/visible absorptivity than their bulk counterparts, they can not be 
the 21 fim feature carrier since the minimum required TiC dust mass is larger than the maximum 
available TiC mass (m™^ ~ 2.5 x 10"^ see §2) by over two orders of magnitude! 

Hony et al. (2003) assumed a hypothetical absorption efficiency for TiC nanograins 



^abs(A) 



AcUV , 912A<A<l^m; 

kUV(2-A) , l^m< A<2/.m; (3) 
, A > 2 fim : 



where is the "constant level" UV/visible mass absorption coefficient. They concluded that 
the astronomical 21 feature could only be accounted for by the TiC model with a reasonable 
amount of TiC dust mass z/nano TiC grains have an UV/visible absorption efficiency as high as 
QUV ^ 8 ^here QUV ^ (4/3) apT^.c^^^. 
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Fig. 2. — The TiC mass required to account for the observed 21 emission feature as a function of the UV/visible 
mass absorption coefficient k^,^ using the Ifony et aL (2003) hypothetical formula (see Eq.[3]) or an optimized Drude 
profile (see §4). We also show the TiC mass derived from models using the optical constants of bulk TiC (filled 
square for grains of a = 10 A). The horizontal lines plot the Kramers-Kronig lower limits to the TiC mass inferred 
from the wavelength-integrated absorption cross sections based on the Hony et al. (2003) approximation (dot-dashed 
line) and the optimized Drude profile (dotted line) (both with F — 1.5; see §4 and Eq.[5]). The shaded area indicates 
reasonable values for the total available TiC dust mass (left axis) and the total envelope mass including both gas and 
dust (right axis). 

In Figure 1 we also plot the Hony et al. (2003) hypothetical mass absorption coefficient with 
'^^hs ~ ^-^ ^ 10^ cm^ g^^ . In comparison with the Kabs('^) values calculated from the optical 
constants of bulk TiC, the Hony et al. (2003) formula greatly enhances the ability of TiC dust 
in absorbing the HD 56126 stellar radiation which peaks at A~ 0.45 ^m. In Figure 2 we show the 
minimum TiC dust mass tti^Iq required to account for the observed 21 /um feature as a function of 
'^ate- '^^^ ^^^^ ^^^^ '^ate ^ ^-^ ^ the minimum required TiC dust mass could 

be smaller than the maximum available TiC mass, implying that the TiC model may be tenable. 

4. Constraints from the Kramers-Kronig Relations 

Let C|!;°|:(A) be the total extinction cross sections of TiC dust at wavelength A, and C*°{:(A)dA 
be the extinction integrated over the entire wavelength range from to oo. As shown by Purcell 
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Fig. 3. — The F(eo; shape) factor as a function of the axial ratio a/b for highly conducting spheroids (i.e., the static 
dielectric constant eo — > oo). For modestly elongated (a/6 < 3) or flattened {b/a < 3) grains, the F factor is not 
expected to significantly deviate from unity. 

(1969), the Kramers-Kronig (KK) dispersion relations can be used to relate /q°° C|!;°j:(A)(iA to the 
total grain volume Vtot through 

/■oo 

/ C*°^(A)dA = 3^VtotF(eo; shape) , (4) 
Jo 

where F is the orientationally-averaged polarizability relative to the polar izability of an equal- 
volume conducting sphere (Purcell 1969; Draine 2003). 

We can apply Eq.(4) to the HD 56126 circumstellar envelope to obtain a lower bound on the 
TiC dust mass m^^, taking the envelope to be a vacuum sparsely populated by spheroidal TiC 
grains.^ Since Cl°l = (C*^* + > C*^* > (where C*^* and C*°* are the total absorption 

and scattering cross sections, respectively), the integration of C*^* over a finite wavelength range 
represents a lower limit to Cl'^{X)dX, and implies a lower limit to the volume of space which 



This is justified since the total volume of the TiC grains 9.0 x 10^^ cm'^ for a total mass of 2.2 x 10~^ JTI©) 
is negligible compared to the total volume of the dusty circumstellar envelope around HD 56126 (~3.1 x 10^^ cm^). 
Apparently, the envelope also satisfies the prerequisite [e — 1| ^ 1 for the KK relations to be applicable (where e is 
not the dielectric constant of TiC, but of the envelope itself). 
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must be filled by TiC grains and hence a lower limit to the total TiC dust mass 

/*oo 

Stt^F (eo; shape) 7912 X ^""'^ ' 

f oo 



TiC 



^TiCfT 

3vr2F ./912A 



/•oo 

— / „ Kabs(A)<iA 
J912A 

(5) 



2pTicg^g^ (21 ;^m) (rg,,, - rg^^J j^^^l ^abs(A)dA 



SttV^ In (rmax/rmin) F J°° ° Kabs (A) F* dX 



•'912A 

The dimensionless factor F depends only upon the grain shape and the static (zero-frequency) 
dielectric constant eo of the grain material (Purcell 1969; Draine 2003). Since TiC is a metallic 
material, we calculate the F factors expected for highly conducting materials (i.e., eg — > oo) as 
a function of the axial ratio a/h (where a and h are the semiaxis along and perpendicular to the 
symmetry axis of the spheroid, respectively). As shown in Figure 3, F would not appreciably 
exceed unity unless the highly conducting grains are highly elongated (for prelates) or flattened 
(for oblates). If they are just modestly elongated ot flattened like interstellar grains (a/6 ~3 
[Greenberg & Li 1996] or h/a ~2 [Lee k, Draine 1985]), we would expect F < 1.5. Therefore, it is 
reasonable to adopt F = 1.5 in the following discussions. 

It is immediately seen in Eq.(5) that one cannot relax the TiC mass requirement just by 
increasing the UV/visible absorption level; instead, the Kramers-Kronig lower limit to the total 
TiC dust mass ?ti^Jq is independent of the absolute level of ^Cabs- As discussed in §3, the nano-TiC 
model appears to be tenable if one adopts the Hony et al. (2003) /tabs (A) formula with > 
5.6 X lO^cm^g"^. However, the KK relations (see Eq.[5]) results in a lower bound of mi^?Q ^ 
2.2 X 10~^ m© to the total TiC dust mass, exceeding the maximum available TiC mass of m™^ ^ 
2.5 X lO~^m0 (see §2)^ by a factor of ~90. Therefore, it appears that the nano-TiC model 
encounters great difficulty in meeting the TiC abundance constraint.^ 

One may argue that the Hony et al. (2003) formula (Eq.[3]) does not make economical use of 
the TiC absorption. A fine tuning on the spectral dependence of Kabs(A) might be able to enhance 
the grain's "ability" to absorb stellar radiation, resulting in an increase in j^^^Ksx^^^X) F^dX 
and a decrease in (see Eq.[5]). We have tried models with a Drude profile- like UV/visible 



^ To be really generous, we can obtain an upper limit of mxic ^ 6.1 x 10"'^ m© by assuming that all Ti elements in 
HD 56126 (with a zcro-agc main-scqucncc mass of ~ 1.1 me) have condensed in the form of TiC. The TiC model with 
Kab^ ^ 2.3 X 10*' cm^ g~^ (see §3 and Eq.[3]) would then appear tenable since mxic ^ 'WtkJ- But the KK lower limit 
on the TiC mass mxic ^ 2.2 x 10~^ tuq is far in excess of ttitk? unless the TiC nanograins are extremely elongated 
{a/h > 25) or flattened (6/o > 80) so that F > 24. 

® In order for the TiC dust mass m!^?^^ derived from the KK relations not to exceed the maximum available TiC 
mass ruTic 5 the TiC nanograins should be extremely elongated {a/b > 44) for prolates or flattened {b/a > 220) for 
oblates so that F > 59. But it is hard to imagine how such extremely-shaped TiC nanograins could form and survive 
in circumstellar environments. 
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absorption spectrum^ i^absW = '^abs (7'^)^/ ('^^ ~ -^o)^ + (7-^)^ with Aq = 0.44 /xm and 7 = 
0.35 /xm (see Fig. 1) ^Hailored^ to optimize the "abihty" of TiC nanograins to absorb the HD 56126 
stellar radiation. But it is found that models using this Kabs functional form require a Kramers- 
Kronig lower limit of « 1.2 x lO~^m0, still ~50 times larger than the maximum available 

TiC mass.^ 



5. Discussion 

In this Letter we have assumed that the TiC dust follows the spatial distribution of the bulk 
carbon dust (see §2). But this is not critical. Even if we assume that all the TiC dust accumulates 
at the inner edge of the envelope, to explain the observed 21 ij,m feature, one still needs ~ 42% of 
the total TiC mass derived in §§3,4, since for the latter the average starlight intensity to which 
the dust is exposed is just ^^2% higher. Therefore, local density enhancement of TiC would not 
overcome the TiC abundance problem. High-resolution (FWHM 0.4") imaging at 20.8 /um revealed 
two bright blobs at ~ 1" from the star, but their surface brightness is at most two times higher 
than the average (Kwok, Volk, h Hrivnak 2002), suggesting the insignificance of local density 
enhancement of the 21 /xm feature carrier. 

Von Helden etal. (2000) postulated that TiC nanograins could have formed during the extreme 
conditions (i.e., high density and high pressure) associated with the short, high mass- loss superwind 
phase where AGB stars lose the remaining stellar envelope, terminating their life on the AGB and 
starting the transition to the PN phase. However, Kwok et al. (2002) argued that there is no 
evidence that the 21 //m emission is created by a sudden ejection at the end of the AGB, implying 
that the necessary physical conditions required to form TiC in the ejecta of carbon-rich evolved 
stars may not be met. 

In summary, although we are lack of experimental knowledge of the UV/visible absorption 
properties of nano TiC, a physical argument based on the Kramers-Kronig dispersion relations 
implies that the nano-TiC model requires too much Ti to reconcile with the observed photospheric 
Ti abundance.^ Recently, the nano-TiC model was also challenged by Chigai et al. (2003), who 
found that the abundance ratio of Ti to Si needed to reproduce the observed flux ratios of the 



'^Spheroidal metallic nanograins are expected to have a Drude-profile like resonance if it is mainly due to free 
electrons (Bohren & Huffman 1983). 

*A narrower Kabs(A) profile leads to a smaller m™, provided its peak is fixed at Ao « 0.44 /um. However, even if 
we adopt an unphysical (5-function for Kabs(A), the corresponding KK lower limit mxix; ~ 6.8 x 10~^ m© still exceeds 
mxic by a factor of ~ 27. We note that the grain surface scattering of electrons alone would result in a plasma 
resonance width of ~ 0.1 /xm for a 10 A grain a typical Fermi velocity vf = 10 cms (see Li 2004). 

®If these grains are coated by a mantle of graphite (see Chigai et al. 2003), the situation would be even worse 
since the large IR emissivity of graphite together with the fact that the increase in grain size may prohibit them from 
single-photon heating would jointly suppress emission in the 21 /um feature. 
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11.3 /Liin SiC feature to the 21 fj,m. TiC feature of post-AGB stars must be at least 5 times larger 
than the solar abundance ratio. Finally, we should note that laboratory absorption spectra of TiC 
nanoclusters also show prominent features at other wavelengths which arc not seen in post-AGB 
stars (e.g. Tii4Ci3 also has a strong band at ~ 16 //m; van Heijnsbergen et al. 1999). 

1 thank D. Arnett, A. Burrows, A.G.W. Cameron, S. Hony, D.S. Lauretta and K. Volk for help- 
ful discussions/comments, T. Koide and D.W. Lynch for providing me with their optical constants 
of bulk TiC. I also thank the University of Arizona for the "Arizona Prize Postdoctoral Fellowship 
in Theoretical Astrophysics" . 
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